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Abstract 

A  method  for  remote  optieal  measurements  of  magnetie  field  variations  above  the  surfaee  of 
seawater  is  analyzed.  This  magnetometry  meehanism  is  based  on  the  faet  that  polarized  laser 
light,  in  the  presenee  of  the  Earth’s  magnetie  field,  will  exhibit  a  polarization  rotation  when 
refieeted  off  the  surfaee  of  the  water  and  off  an  underwater  objeet.  The  two  meehanisms 
responsible  for  the  polarization  rotation  are  the  Surfaee  Magneto-Optieal  Kerr  Effeet  (SMOKE) 
and  the  Earaday  effeet.  In  both  meehanisms,  the  degree  of  polarization  rotation  is  proportional  to 
Earth’s  loeal  magnetie  field.  Variations  in  the  Earth’s  magnetie  field  due  to  an  underwater  objeet 
will  result  in  variations  in  the  polarization  rotation  of  the  laser  light  refieeted  off  the  water’s 
surfaee  (SMOKE)  and  off  the  underwater  objeet  itself  (Earaday  effeet).  An  analytieal  expression 
is  obtained  for  the  polarization-rotated  field  when  the  ineident  plane  wave  is  at  arbitrary  angle 
and  polarization  with  respeet  to  the  water’s  surfaee. 
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I.  Introduction 

Optical  magnetometry  can  be  a  highly  sensitive  method  for  measuring  small 
variations  in  magnetic  fields  [1-3].  The  development  of  a  remote  optical  magnetometry  system 
would  have  important  applications  for  the  detection  of  underwater  and  underground  objects  that 
perturb  the  loeal  ambient  magnetic  field.  For  a  number  of  magnetie  anomaly  detection  (MAD) 
applications,  juG  magnetic  field  variations  must  be  detected  at  standoff  distances  greater  than 

lO^mfrom  the  sensor  [4-5].  The  remote  atmospherie  optieal  magnetometry  mechanism 
considered  here  is  based  on  polarization  changes  in  reflected  laser  light  from  the  surface  of  sea 
water. 

The  dielectrie  properties  of  water  in  the  presence  of  a  magnetic  field  will  rotate  the 
polarization  of  an  optical  field.  The  polarization  rotation  oceurs  via  two  distinct  mechanisms: 
the  Surface  Magneto-Optical  Kerr  Effect  (SMOKE)  and  the  Earaday  rotation  effect.  SMOKE  is 
purely  a  surface  phenomenon  in  which  the  degree  of  polarization  rotation  of  the  reflected  light  is 
proportional  to  the  magnetic  field  but  not  to  the  propagation  distance  in  the  water  [6].  The 
SMOKE  mechanism  is  used  in  material  science  researeh  in  such  devices  as  the  Kerr  microscope. 
The  Earaday  effect  is  due  to  the  difference  in  propagation  speed  of  right-hand  and  left-hand 
polarized  light,  and  the  resulting  degree  of  polarization  rotation  is  proportional  to  both,  the 
magnetic  field  and  the  propagation  distanee. 

In  this  paper,  we  consider  whether  these  effects  can  be  used  for  the  remote  deteetion 
of  underwater  objeets.  An  underwater  object  can  produce  appreeiable  magnetic  field 
perturbations  above  the  water’s  surface.  The  proposed  magnetometry  meehanism  is  shown 
sehematically  in  Eig.  1.  Here,  a  linearly-polarized  laser  beam  is  propagated  to  and  reflects  off 
the  water’s  surface  near  the  detection  site.  The  degree  of  polarization  rotation  in  the  reflected 
field  is  proportional  to  the  loeal  magnetic  field.  The  magnitude  of  the  polarization-rotated 
eomponent  of  the  reflected  field  is  measured  at  various  loeations  near  the  inspection  site.  An 
irregularity  in  the  measurement  at  one  loeation  would  eorrespond  to  a  deviation  in  the  local 
magnetic  field.  The  local  field  deviation  eould  indicate  the  presence  of  an  underwater  objeet. 
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Figure  1 :  Polarized  laser  light  is  propagated  to  the  surfaee  of  the  water  and  exhibits  polarization 
rotation  upon  reflection  due  to  Faraday  rotation  and  SMOKE.  The  purpose  of  this  paper  is  to 
estimate  the  polarization  rotation  angle  A0  cc  .  If  A0  is  large  enough  to  be  precisely 
measured,  the  mechanism  shown  in  this  figure  enables  the  measurement  of  fluctuations  in  the 
Earth’s  magnetic  field  ,  which  in  turn  may  indicate  the  presence  of  an  underwater  object. 
The  solid  arrows  in  the  figure  denote  wave  vectors. 


II.  Model 

To  analyze  the  Surface  Magneto-Optical  Kerr  and  the  Faraday  rotation  effects,  we 
first  obtain  the  dielectric  properties  of  magnetized,  pure  water.  The  polarization  field  associated 
with  magnetized  water  at  the  laser  frequency  is  obtained  using  an  extended  form  of  the  Lorentz 
model  [7].  The  polarization  field  is  used  to  derive  a  wave  equation  for  the  reflected  laser  field. 
By  appropriately  matching  boundary  conditions  at  the  various  interfaces,  we  obtain  the  reflected 
field  which  has  its  polarization  modified  by  both,  the  SMOKE  and  Faraday  effect.  The 
polarization-rotated  component  of  the  field  is  proportional  to  the  local  magnetic  field.  In 
applying  the  boundary  conditions,  the  laser  propagation  direction  and  polarization  as  well  as  the 
Earth’s  magnetic  field  are  taken  to  be  in  an  arbitrary  direction. 
i)  Polarization  Field  of  Magnetized  Water 

In  the  Eorentz  model,  the  dielectric  properties  of  the  molecule  are  represented  by  the 
motion  of  the  forced  displacement  of  the  electron  distribution  from  its  equilibrium  position.  The 
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displacement  of  the  eleetron  distribution  5r{r,t)  is  modeled  by  the  following  driven  harmonie 
oseillator  equation, 


— j  +  y —  5r  =  —  E+ - xB„/c  , 

dt^  dt]  m  dt  " 


where  is  the  electron  binding  frequency,  y  is  the  damping  coeffieient,  q  is  the  eleetronic 
charge,  m  is  the  electronie  mass,  andfi^  is  the  loeal  magnetie  field,  i.e..  Earth’s  magnetie  field 

-0.5  G.  The  laser  eleetric  field  is  E(r,t)  =  Re(^E(r)exp(-/(yt)^  where  u;  is  the  frequeney. 

The  magnetie  field  assoeiated  with  the  laser  is  negleeted  on  the  right  hand  side  of  Eq.  (1)  since  it 
produees  forees  at  frequencies  of  Ico  and  0. 

The  polarization  field  assoeiated  with  the  magnetized  water  moleeules  is  given  by 

P  (r)  =  qNdf  ,  where  Sfir)  represents  the  eleetronic  displacement  from  equilibrium  of  the 

water  moleeule,  N  is  the  effective  bound  eleetron  density  of  the  water  molecules,  and  phasor 
notation  is  used  throughout  to  represent  the  various  field  quantities,  i.e., 

e(r,t)  =  Re(0(r)  exp(-/(ut)j .  The  displacement  5f{r)  is  given  by 

-  iyct>jSf(r)  =  q(^E(r)  -  ict>Sf(r)xB^  /  c^/ m.  (2) 

The  polarization  field,  eorreet  to  first  order  in  the  magnetic  field,  is 

2 

P(r)  =  qNSf^  xE(r)/D(<o)),  (3) 

where Z)((y)  =  -  0/  -  iyco,  =  [Anq^N / ,  and  Qq  =  qB^  / mc~  10^5  ' . 

The  laser  frequency  is  chosen  to  be  in  the  blue-green  range  where  the  penetration  depth  is 
maximum.  In  the  following  examples,  the  laser  frequeney  is  taken  to  be  ®  =  3.8x10'^ rad/s 
(ru  =  2.48eV /  ti)  where  pure  water’s  absorption  ooeffieient  is  minimum,  i.e.,  a  =  5 x lO^cm"' ; 
the  eorresponding  penetration  depth  (e-folding  length)  is  «  20  m  .  In  the  blue-green  regime, 

the  real  part  of  the  refractive  index  is  Re[Uo] » 1.34 .  The  lowest  energy  electron-binding 
frequency  in  water  is  Qg  ~  23.5eV /  h ,  »  21eV /  h ,  and  the  damping  coefficient  is 

y  «  3x10‘^5  ' . 
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In  general,  there  are  many  resonanees  above  whieh  eontribute  to  the  refraetive  index  of 
water.  We  negleet  these  binding  frequeneies  sinee  they  are  at  higher  energies.  In  addition,  the 
oseillator  strength  parameter  assoeiated  with  the  Qg  resonanee  is  taken  to  be  unity  [8]. 

Nonlinear  effeets  in  water  ean  be  negleeted  sinee  the  laser  intensity  is  well  below  the  level  where 
these  beeome  important. 
ii)  Wave  Equation 

In  the  absenee  of  free  eharges  and  eurrents.  Maxwell’s  equations  for  the  fields  in 
phasor  notation  are  given  by  VxH  =  -i(ft»/c)D,  V-D  =  0 ,  VxE  =  i((y/c)B  andV-B  =  0. 
These  equations  ean  be  eombined  to  give 

V'E  +  (0"/c")D  =  V(V-E).  (4) 


To  obtain  a  wave  equation  for  E  ,  we  note  that  the  eleetrie  flux  density  D  =  (E  +  4;r  P)  is  given 
by  D  =  +  i{s^-\){a)l  D{a)))il^  x  )E,  where  the  eleetrie  polarization  field  is  given  in  Eq.  (3) 

and  where  the  dieleetrie  eonstant  of  water,  in  the  absenee  of  the  Earth’s  magnetie  field,  is 
=  1  +  /  D{co) .  Sinee  V  •  D  =  0 ,  we  find  that 

V .  E  =  f  V-(n.  X  E)  =  “  (il,  ■  (V  X  E)) ,  (5) 

D{a))  D(co) 

where  the  final  form  for  V  •  E  is  valid  sinee  is  spatially  uniform.  Using  Eqs.  (4)  and  (5),  the 
wave  equation  takes  the  form 


we  +  ^„4e  =  -^-4(^.-i) 


,ry  ®  ^  V(a^-(VxE))^ 


c"  c~  D(ci))  D(co) 

The  right-hand  side  of  Eq.  (6)  ean  be  simplified  by  noting  that 

•  (V  X  E))  =  X  (V  X  V  X  E))  -i-  •  V)(V  x  E)  and,  to  lowest  order  in  , 

V  X  V  X  E  =  .  The  final  form  for  the  wave  equation,  eorreet  to  first  order  in  ,  is 


(6) 


,2  A 


^  +^o  — 

V  C  y 


E  =  ^  (a„-V)(VxE)  . 


(7) 


7)(co) 

Sinee  the  right-hand  side  of  Eq.  (7)  is  small,  we  ean  write  the  total  field  as  E(r)  =  E^  -l-  (7E , 
where  E^  here  denotes  the  laser  field  in  the  absenee  of  the  perturbation  (i.e.  the  Earth’s  magnetie 
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field)  and  is  the  small  polarization-modified  eomponent  of  the  laser  field  due  to  the  Earth’s 
magnetie  field  and  driven  by  the  laser  field  .  Separate  wave  equations  ean  be  written  for 
these  fields: 

V"  +  «Ae„=0,  (8a) 

I  c  J 

+  g,  4 W  =  f ~  ‘ f  (a.  ■  V)( V  X  E. ) .  (8b) 

I  C  )  D((o) 

In  the  special  case  where  the  laser  field  E^  is  transverse  to  the  magnetic  field  ,  one  can 
recover  from  Eq.  (8b)  the  usual  Earaday  rotation  angle,  i.e. 

A6p  w  o/Q.^LI (2c  Qg)  □  (lO^G^^m  4ad)L5^ ,  where  L  is  the  propagation  distance  [9]. 

These  wave  equations  are  numerically  solved  for  arbitrary  laser  field  propagation  direction, 
polarization  and  magnetic  field  direction. 

III.  Analysis  of  Polarization  Rotation 

The  proposed  magnetometry  method  is  feasible  only  if  the  laser  beam’s  polarization 
rotates  a  substantial  amount  upon  reflection  off  the  water’s  surface.  The  purpose  of  the  analysis 
is  thus  to  calculate  the  amount  of  polarization  rotation  A0  =  |^£'|/|£’^j|,  where  AE  is  the 

polarization-rotated  component  of  the  reflected  laser  field  and  is  the  un-rotated  component. 
Alternatively,  this  is  written  as  A0  =  (l/R)|AE|/|Eg|  ,  where  Eg  is  the  incident  laser  field  and  the 
reflection  coefficient  is  R  =  (n  -  l)/(n  -i- 1)  »  0. 1 5 . 

We  first  consider  a  simplified  configuration  as  shown  in  Eig.  2  to  obtain  estimates 
for  the  amount  of  polarization  rotation  A0 .  Here  a  linearly-polarized,  monochromatic  laser 
beam  is  propagated  parallel  to  Earth’s  magnetic  field  which  is  taken  to  be  perpendicular  to  the 
surface  of  the  water.  An  optically  reflective  underwater  object  at  depth  h  below  the  water’s 
surface  is  represented  by  a  perfectly  conducting  plate  parallel  to  the  surface. 

Erom  this  simplified  configuration  will  emerge  two  distinct  contributions  to  the 
polarization  rotation:  the  Earaday  effect  and  SMOKE.  The  Earaday  effect  yields  the  ratio 
\5EjE^~\2nnQh^l  X ,  whereas  SMOKE  yields  |AE/Eg|  3x10''^ .  The  Earaday  effect  is 
larger  by  the  ratio  h  /  /I  □  1  . 
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The  SMOKE  contribution  to  the  polarization  rotation  is  obtained  for  an  arbitrary 
configuration. 
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Figure  2:  A  simple,  geometrical  configuration  of  reflected  laser  light  undergoing  polarization 
rotation  due  to  Faraday  rotation  and  SMOKE.  The  purpose  of  this  configuration  is  to  obtain 
estimates  for  the  polarization  rotation  angle  IS.6 . 


The  solid  arrows  in  Fig.  2  are  the  wave  vectors  while  the  dashed  arrows  denote  the  electric  fields 
corresponding  to  the  incident  wave  E^,  the  transmitted  wave  Ej. ,  the  object-reflected  wave  E^2  ’ 
and  the  surface-reflected  wave  E^j  which  is  modified  by  a  small  polarization-rotated  component 
SE .  The  electric  field  can  be  written  as, 

E(z)  =  Eg  exTp(ikz)e^  +  exp(-ikz)e^ +SE  exp(-ikz)e^  ,  ( z  <  0)  (9a) 

E (z)  =  Ej  exp(ik^z)e^  +  Ej.  exp(ik_z)e_  + 

exp(-i^_^z)e^  -h  E^2  exp(-i^_z)e_ ,  (0<z<h) 

where  =(e^±ie^.)  /  yfl  ,  and  are  unit  vectors  in  the  x  and  y  direction  respectively, 
k  =  0)/  cis  the  wavenumber  in  air,  =ck^  /  is  the  complex  index  of  refraction  in 

magnetized  water,  and  =(£„  -  1)®  Q/  Z)(®)  .  Application  of  the  electromagnetic  boundary 
conditions  yields  for  the  configuration  shown  in  Fig.  2  the  following  ratio  for  the  polarization- 
rotated  field. 
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In  the  limit  that  /z  — >  oo ,  the  re  fleeted  Faraday  rotation  effect  vanishes  and  the  polarization 
rotation  is  due  solely  to  SMOKE.  In  this  limit  Eq.  (10)  reduces  to 


5E 


—n_ 

(l  +  nj(l  +  n_) 


(11) 


where  »(no-l)(ft)Qo/Qg)/(no +no)  □  3x10  .  The  ratio  of  the  fields  given  in  Eq.  (10)  is 

shown  in  Eig.  3  as  a  funetion  of  the  depth  of  the  eonducting  plate  h  . 


Eigure  3:  Ratio  of  the  polarization-rotated  refleeted  fleld  to  the  ineident  field  as  a  funetion  of 
object  depth  h  .  The  polarization  rotation  results  from  a  combination  of  Earaday  rotation  and 
SMOKE;  however,  the  eontribution  due  to  SMOKE  is  small. 


We  now  consider  the  general  SMOKE  contribution  to  the  polarization  rotation  by 
analyzing  the  proeess  for  arbitrary  field  directions.  The  polarization-rotated  field  due  to  SMOKE 
as  a  funetion  of  the  field  eomponents  of  an  ineident  plane  wave,  for  arbitrary  polarizations  and 
angle  of  ineidence  as  well  as  arbitrary  orientation  of  the  magnetic  field,  is  obtained  in  Appendix 
A  and  found  to  be, 

\5E\  □  K(6») \n,  I  (eIe]^  +  2  [eIe]^  +  )  cos'  6 + EJe;'  cos'  o")"  ,  (13) 
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where  K(0)  ^  (n^(i>/mnl)g(0)  =  2.5x10 □  2.5x10  and 

rig  w  (cosdg  -  ^sin^  Og  sin^  / 1.8)  ~  1  .  In  Eq.  (12),  6  is  the  angle  of  the  incident  wave  vector 
with  respect  to  the  water’s  surface  normal  vector.  The  angles  6g  and  (j)g  are,  respectively,  the 
polar  and  azimuthal  angle  of  Earth’s  magnetic  field  with  respect  to  the  incident  wave  vector. 
Appendix  A. 


m 


Eigure  4;  The  function  g{0)  plotted  above  represents  the  angular  dependence  of  the 
polarization-rotated  field  on  the  incident  wave  vector’s  orientation  for  a  specific  geometry 
{9g=(l>g=Q  ,  E^  =Q). 

We  find  that  no  choice  of  polarization  or  angle  of  incidence  will  significantly 
increase  the  SMOKE  contribution  to  the  polarization  rotation:  see  Eig.  4. 

IV.  Discussion 

A  mechanism  for  the  remote  optical  detection  of  a  shallow,  i.e.,  -  1  -  10  m, 
underwater  conducting  object  has  been  presented.  In  this  mechanism,  a  polarized  blue-green 
laser  beam  is  reflected  off  the  water’s  surface  and  the  reflected,  polarization-rotated  electric  field 
recorded.  A  variation  in  the  measurements  at  two  nearby  locations  may  indicate  the  presence  of 
an  underwater  conducting  object.  The  polarization  rotation  is  mainly  due  to  the  Earaday  effect; 
SMOKE  has  a  negligible  contribution. 

As  an  alternative  to  the  mechanism  proposed,  one  can  detect  the  underwater  object 
directly  by  measuring  the  un-rotated  component  of  the  reflected  laser  field.  This  reflected 
component  can  be  significantly  larger  in  amplitude  than  the  polarization-rotated  component.  By 
time  gating  the  laser  pulses,  the  object’s  depth  can  be  determined. 
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Appendix  A 

A  perturbative  solution  is  deseribed  here  for  the  polarization-rotated  reflected  field 
associated  with  SMOKE.  The  treatment  allows  for  the  incident  plane  wave  to  have  arbitrary 
angle  and  polarization  with  respect  to  the  water’s  surface  and  allows  for  arbitrary  magnetic  field 
orientation.  The  relevant  coordinate  systems  are  labeled  in  Fig.  Al. 


Figure  Al;  Coordinate  systems  used  in  the  calculation  of  the  polarization-rotated  reflected  field 
associated  with  SMOKE  given  an  arbitrary  geometric  configuration. 


The  arrows  correspond  to  the  transmitted  wave  vector  in  the  absence  of  a  local  magnetic  field 
and  the  incident  wave  vector;  the  boundary  conditions  ensure  that  they  are  coplanar  and  that 
sin  a  =  ((Z»/c)sin^  ,  where  is  the  magnitude  of  the  transmitted  wave  vector.  The 
coordinate  systems  are  defined  such  that  x  =  x'-x”  coincide. 

The  fields  E{y)  ,  D{y)  ,  B(r)=-i(c  /  ry)V  xE(r)  are  expanded  to  first  order  in  the 
small  parameter  A  =  cl)Q,q  /  T>(®)  □  10  "  where,  as  before,  D[co)  =  -  a/  -  iyco  .  Since 

V.D(r)=0  ,  to  each  order  in  A  ,  Eq.  (4)  yields  the  wave  equation  /  c^)Eq  =  0 , 

where  the  field  subscript  denotes  zeroth  order  in  A  .  This  wave  equation  is  used  to  calculate  the 
zeroth  order  field  amplitudes. 

For  z '  >  0 ,  Eq.  (7)  implies  that,  for  the  transmitted  field  Ej.{r)=Ej.  o(r)  -1-  SEj.{r) , 
+co^  Sq!  c^'^dE^{r)  =  i{^£^-\)co{p.^y^V^{v  y.E^  g(r)^/(£:o  T)(®))  (13) 

The  field  is  rewritten  as  (5'E^(r)=r(r)exp(i/:oz")=r(r)epexp(i/:oz")  in  order  to  employ  the 
paraxial  approximation  -1- lik^d /  dz” ,  which  states  that  the  field  varies  a  negligible 
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amount  for  transverse  displaeements  Arj^  on  the  order  of  a  wavelength.  The  boundary 
eonditions  dictate  that  r(r)hasno  x'or  y ' dependence  at  z'  =  0,i.e.  r(r)=r(z').  Applying  the 
paraxial  approximation  to  Eq.  (13)  and  integrating  over  the  coordinate  z'  ,  one  obtains 
r(z')  =  a+bz'  .  The  constant  of  integration  a  is  determined  by  the  continuity  of  the  magnetic 

field  B  =  -/  (c  /  ry)  (V  X  E)  at  the  boundary.  Simplification  of  the  expression  for  yields 

SEj(z'  =  0)  ,  i.e.  \SE\  as  given  by  Eq.  (12). 
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